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We consider the entanglement dynamics between two Unruh-DeWitt detectors at rest separated 
at a distance d. This simple model when analyzed properly in quantum field theory shows many 
interesting facets and helps to dispel some misunderstandings of entanglement dynamics. We find 
that there is spatial dependence of quantum entanglement in the stable regime due to the phase 
difference of vacuum fluctuations the two detectors experience, together with the interference of 
the mutual influences from the backreaction of one detector on the other. When two initially 
entangled detectors are still outside each other's light cone, the entanglement oscillates in time with 
an amplitude dependent on spatial separation d. When the two detectors begin to have causal 
contact, an interference pattern of the relative degree of entanglement (compared to those at spatial 
infinity) develops a parametric dependence on d. The detectors separated at those d with a stronger 
relative degree of entanglement enjoy longer disentanglement times. In the cases with weak coupling 
and large separation, the detectors always disentangle at late times. For sufficiently small d, the 
two detectors can have residual entanglement even if they initially were in a separable state, while 
for d a little larger, there could be transient entanglement created by mutual influences. However, 
we see no evidence of entanglement creation outside the light cone for initially separable states. 

PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.-a 



I. INTRODUCTION 



Recently we have studied the disentanglement process between two spatially separated Unruh-DeWitt (UD) de- 
tectors (pointlike objects with internal degrees of freedom) or atoms, described by harmonic oscillators, moving in a 
common quantum field: One at rest (Alice), the other uniformly accelerating (Rob) These two detectors are set to 
be entangled initially, while the initial state of the field is the Minkowski vacuum. In all cases studied in [l[ , we obtain 
finite-time disentanglement (called "sudden death" of quantum entanglement 0]), which are coordinate dependent 
while the entanglement between the two detectors at two spacetime points is independent of the choice of time slice 
connecting these two events. Around the moment of complete disentanglement there may be some short-time revival 
of entanglement within a few periods of oscillations intrinsic to the detectors. In the strong-coupling regime, the 
strong impact of vacuum fluctuations experienced locally by each detector destroys their entanglement right after the 
coupling is switched on. 

In the above situation we find in [l( the event horizon for the uniformly accelerated detector (Rob) cuts off the 
higher-order corrections of mutual influences, and the asymmetric motions of Alice and Rob obscure the dependence 
of the entanglement on the spatial separation between them. To understand better how entanglement dynamics 
depends on the spatial separation between two quantum objects, in this paper we consider the entanglement between 
two detectors at rest separated at a distance d, possibly the simplest setup one could imagine. This will serve as a 
concrete model for us to investigate and explicate many subtle points and some essential misconceptions related to 
quantum entanglement elicited by the classic paper of Einstein-Podolsky- Rosen (EPR) Q . 
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A. Entanglement at spacelike separation: quantum nonlocality? 

One such misconception (or misnomer, for those who understand the physics but connive to the use of the ter- 
minology) is "quantum nonlocality" used broadly and often too loosely in certain communities [23l |. Some authors 
think that quantum entanglement entails some kind of "spooky action at a distance" between two spacelike separated 
quantum entities (qubits, for example), and may even extrapolate this to mean "quantum nonlocality." The phrase 
"spooky action at a distance" when traced to the source 0] refers to the dependence of "what really exists at one 
event" on what kind of measurement is carried out at the other, namely, the consequence of measuring one part of an 
entangled pair. Without bringing in quantum measurement, one cannot explore fully the existence or consequences of 
"spooky action at a distance" but one could still talk about quantum entanglement between two spacelike separated 
qubits or detectors. This is the main theme of our present investigation. We show in a simple and generic model with 
calculations based on quantum field theory (QFT) that nontrivial dynamics of entanglement outside the light cone 
does exist. 

Another misconception is that entanglement set up between two localized quantum entities is independent of their 
spatial separation. This is false for open systems interacting with an environment (24[. This has already been shown 
in two earlier investigations of the authors [3, Q and will be again in this paper. 

A remark on nonlocality, or lack thereof, in QFT is in place here. QFT is often regarded as "local" in the sense that 
interactions of the fields take place at the same spacetime point [25|, e.g., for a bosonic field <f>(x), a local theory has no 
coupling of 4>{x) and <fr(y) at different spacetime points x and y. It follows that the vacuum expectation value of the 
commutator ( [4>(x), 4>{y)] ) vanishes for all y outside the light cone of x, which is what causality entails. Nevertheless, 
the Hadamard function ( {<fi(x), 4>{y)} ) is nonvanishing in general, no matter x — y is spacelike or timelike. In physical 
terms the Hadamard function can be related to quantum noise in a stochastic treatment of QFT [9( . In this restricted 
sense one could say that QFT has certain nonlocal features. Of course it is well known that in QFT processes occurring 
at spacelike separated events such as virtual particle exchange are allowed. 

B. Issues addressed here 

With a careful and thorough analysis of this problem we are able to address the following issues: 

1) Spatial separation between two detectors. -Ficek and Tanas [l2| as well as Anastopoulos, Shresta, and Hu (ASH) 
Q studied the problem of two spatially separated qubits interacting with a common electromagnetic field. The former 
authors while invoking the Born and Markov approximations find the appearance of dark periods and revivals. ASH 
treat the non-Markovian behavior without these approximations and find a different behavior at short distances. In 
particular, for weak coupling, they obtain analytic expressions for the dynamics of entanglement at a range of spatial 
separation between the two qubits, which cannot be obtained when the Born-Markov approximation is imposed. A 
model with two detectors at rest in a quantum field at finite temperature in (l+l)-dimensional spacetime has been 
considered by Shiokawa in fl3j , where some dependence of the early-time entanglement dynamics on spatial separation 
can also be observed. 

In [l[ we did not see any simple proportionality between the initial separation of Alice and Rob's detectors and the 
degree of entanglement: The larger the separation, the weaker the entanglement at some moments, but stronger at 
others. We wonder if this unclear pattern arises because the spatial separation of the two detectors in [l[ changes in 
time and also in coordinate. In our present problem the spatial separation between the two detectors is well defined 
and remains constant in Minkowski time, so the dependence of entanglement on the spatial separation should be much 
clearer and distinctly identifiable. 

2) Stronger mutual influences. -Among the cases we considered in the largest correction from the mutual in- 
fluences is still under 2% of the total while we have only the first and the second-order corrections from the mutual 
influences. There the difficulty for making progress is due to the complicated multidimensional integrations in com- 
puting the back-and-forth propagations of the backreactions sourced from the two detectors moving in different ways. 
Here, for the case with both detectors at rest, the integration is simpler and in some regimes we can include stronger 
and more higher-order corrections of the mutual influences on the evolution of quantum entanglement. 

3) Creation of entanglement and residual entanglement.— In addition to finite-time disentanglement and the revival 
of quantum entanglement for two detectors initially entangled, which have been observed in [l| for a particular initial 
state, we expect to see other kinds of entanglement dynamics with various initial states and how it varies with spatial 
separations. Amongst the most interesting behavior we found the creation of entanglement from an initially separated 
state [14j and the persistence of residual entanglement at late times for two close-by detectors |l5[ . 
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C. Summary of our findings 

When the mutual influences are sufficiently strong (under strong coupling or small separation), the fluctuations of 
the detectors with low natural frequency will accumulate, then get unstable and blow up. As the separation approaches 
a merge distance (quantified later), only for detectors with high enough natural frequencies will the fluctuations not 
diverge eventually but acting more and more like those in the two harmonic oscillator (2HO) quantum Brownian 
motion (QBM) models [TBI, [Tg] (where the two HOs occupy the same spatial location) with renormalized frequencies. 

If the duration of interaction is so short that each detector is still outside the light cone of the other detector, 
namely, before the first mutual influence reaches one another, the entanglement oscillates in time with an amplitude 
dependent on spatial separation: At some moments the larger the separation the weaker the entanglement, but at other 
moments, the stronger the entanglement. While such a behavior is affected by correlations of vacuum fluctuations 
locally experienced by the two detectors without causal contact, there is no evidence for entanglement generation 
outside the light cone suggested by Franson in Ref. [TTf| . 

For an initially entangled pair of detectors, when one gets inside the light cone of the other, certain interference 
patterns develop: At distances where the interference is constructive the disentanglement times are longer than those 
at other distances. This behavior is more distinct when the mutual influences are negligible. For the detectors 
separable initially, entanglement can be generated by mutual influences if they are put close enough to each other. 

At late times, under proper conditions, the detectors will be entangled if the separation is sufficiently small, and 
separable if the separation is greater than a specific finite distance. The late-time behavior of the detectors is governed 
by vacuum fluctuations of the field and independent of the initial state of the detectors. 

Since the vacuum can be seen as the simplest medium that the two detectors immersed in, we expect that the 
intuitions acquired here will be useful in understanding quantum entanglement in atomic and condensed matter 
systems (upon replacing the field in vacuum by those in the medium). To this extent our results indicate that the 
dependence of quantum entanglement on spatial separation of qubits could enter in quantum gate operations (see Q 
for comments on possible experimental tests of this effect in cavity ions), circuit layout, as well as having an effect on 
cluster states instrumental to measurement-based quantum computing. 



D. Outline of this paper 



This paper is organized as follows. In Sec. II we describe our model and the setup. In Sec. Ill the evolution 
of the operators is calculated, then the instability for detectors with low natural frequency is described in Sec. IV. 
We derive the zeroth-order results in Sec. V, and the late-time results in Sec. VI. Examples with different spatial 
separations of detectors in the weak-coupling limit are given in Sec. VII. We conclude with some discussions in Sec. 
VIII. A late-time analysis on the mode functions is performed in Appendix A, while an early-time analysis of the 
entanglement dynamics in the weak-coupling limit is given in Appendix B. 



II. THE MODEL 



Let us consider the Unruh-DeWitt detector theory in (3+l)-dimensional Minkowski space described by the action 
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E 

j=A,B 



drj- +A d A x<S>{x) dTjQjfa)? (x* - z?(tj)) ,(1) 



where the scalar field $ is assumed to be massless, and Ao is the coupling constant. Qa and Qb are the internal degrees 
of freedom of the two detectors, assumed to be two identical harmonic oscillators with mass tuq = 1, bare natural 
frequency Oo, and the same local time resolution so their cutoffs in two-point functions [171 ] are the same. The left 
detector is at rest along the world line z^(t) = (t, —d/2, 0, 0) and the right detector is sitting along z^(t) = (i, d/2, 0, 0). 
The proper times for Qa and Qb are both the Minkowski time, namely, ta = tb = t. 

We assume at t — the initial state of the combined system is a direct product of the Minkowski vacuum | 0m ) 
for the field $ and a quantum state | Qa, Qb ) f° r the detectors Qa and Qb, taken to be a squeezed Gaussian state 
with minimal uncertainty, represented by the Wigner function of the form 



p(Qa, Pa, Qb,Pb) = —ttz cxp -~ 



2 9 1 9 Of 2 9 1 9 

72 (Qa + Qb) +^(Qa- Qb) + 72 {Pa - Pb) +i*>(Pa + Pb) 



(2) 
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How the two detectors are initially entangled is determined by properly choosing the parameters a and (3 in Qa and 
Qb- When (3 2 = h 2 /a 2 , the Wigner function ([5]) becomes a product of the Wigner functions for Qa, Pa and for Qb, 
Pb, thus separable. If one further chooses a 2 = h/tt, then the Wigner function will be initially in the ground state 
of the two free detectors. 

After t = the coupling with the field is turned on and the detectors begin to interact with each other through the 
field while the reduced density matrix for the two detectors becomes a mixed state. The linearity of (pj guarantees 
that the quantum state of the detectors is always Gaussian. Thus the dynamics of quantum entanglement can be 
studied by examining the behavior of the quantity £ jl| and the logarithmic negativity £V [HI : 



£ 
Em 



det 



V PT + —M 



ih. 



max {0, — log 2 2c_} . 



(3) 
(4) 



Here M is the symplectic matrix 1 <8> (—i)a y , V PT is the partial transpose ((Qa, Pa,Qb,Pb) — ► (Qa, Pa,Qb, —Pb)) 
of the covariance matrix 



vaa v A b 
vba v B b 



(5) 



in which the elements of the 2x2 matrices Vy are symmetrized two-point correlators Vy m " = ( 7?.™, 7Z™ ) = ( {TV^TV-- 
1ZJ1IY 1 ) > /2 with Kf = (Q i (t),P l (t)), m,n= 1,2 and i,j = A, B. (c+,c_) is the symplectic spectrum of V PT - 
(iTi/2)M, given by 



Z± ^Z 2 -4detV 



1 1/2 



with 



Z = detvAyi + detv_BB — 2detvAs- 



(6) 



(7) 



For the detectors in Gaussian state, Ej^ > 0, S < 0, and c_ < Ti/2, if and only if the quantum state of the detectors 
is entangled [l9l ]. Ej,/ is an entanglement monotone (2p| whose value can indicate the degree of entanglement: below 
we say the two detectors have a stronger entanglement if the associated E^r is greater. In the cases considered in 
Ref. ml and this paper, the behavior of £ is similar to — E/f when it is nonzero. Indeed, the quantity S can also be 
written as 



detV 
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(8) 



We found it is more convenient to use £ in calculating the disentanglement time. We also define the uncertainty 
function 



T = det 



V + i-M 



(9) 



so that T > is the uncertainty relation [l9( . 

To obtain these quantities, we have to know the correlators ( H™ ,71™ ), so we are calculating the evolution of 
operators 7?.™ in the following. 



III. EVOLUTION OF OPERATORS 

Since the combined system ((T|) is linear, in the Heisenberg picture [UH]]], the operators evolve as 

T 



Qi(t) = 



gi +) (t,k)S k + 5 r 3 (t,k)i>i 



f (+ \x,k)b k +f(-\x,k)hi 



(10) 
(11) 
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with i,j = A,B. qf \ Qi , and are the (c-number) mode functions, dj and are the lowering and 

raising operators for the free detector j, while &k and 6^ are the annihilation and creation operators for the free field. 
The conjugate momenta are Pj(t) = d t Q j (t)a,iid tl(x) — d t Q{x). The evolution equations for the mode functions 
have been given in Eqs.(9) — (12) in Ref. [lj with ZAit) and zb(t) there replaced by z^(t) = (i, — d/2, 0, 0) and 
Zg(t) = (i, d/2,0, 0) here. Since we have assumed that the two detectors have the same frequency cutoffs in their 
local frames, one can do the same renormalization on frequency and obtain their effective equations of motion under 
the influence of the quantum field : 

(d? + 2 1 d t + nl)q ( i >\t) = ^9{t-d)^\t-d), (12) 

(d 2 + 2 1 d t + n 2 r )q^\t 7 k) = ^0(t-<Q^ + >(t-d,k)+Ao/ o (+ W)>k)> (13) 

where q~B = <1A, q_A = <Ib, is the renormalized frequency obtained by absorbing the singular behavior of the 
retarded solutions for and around their sources (for details, see Sec.IIA in Ref. 17]). Also 7 = Aq/87t, and 
/q (a;,k) = e - Mt + ikx ; with uj = |k|. Here one can see that qs and qA are affecting, and being affected by, each 
other causally with a retardation time d. 



luie 



lk-x 



The solutions for q\ and q\ satisfying the initial conditions /( + )(0,x;k) = e lk x , dtf^ (0, x; k) = 
qW(0) = 1, d t qf(0) = -iO r , and /^(0,x) = d t tf\o,x) = g<+>(0;k) = d t q^(0;k) = $\o) = d t qf(0) = (no 
summation over j) are 

qf\k;t) = ^^^^-nd) ^y i e (-i)"*4{(M 1 -M 2 )"+ 1 e-^( t -" d ) + 

71=0 ^ ' 

e _ 7 (t-„ d ) y (Mi _ Ma)n -m [M 2 W m (t - nd) - AW* (t - nd)] 1 , (14) 



m=Q 



and 



= ?f = E^ +1 (15) 

n=0 n=0 

(no summation over j), where = ^/ Sl^ — Mi = (—lu — ij + fi) -1 , M 2 = (—to — 17 — O) -1 , Wo(t) = e int , 

W n (t) = / dt n ^i smf2(t - t n -i) / dt n -2sinfl(t n -i-t n -2)--- dt sin Q(tx - t )W (t ), (16) 
Jo Jo Jo 



for n > 1, and 



q n (t) = 8(t - nd) e -7(*-nc«) [siWn(i _ nd) + S2 ^ (t _ w)] 



(17) 



with si = [1 - n-^O,. + ij)}/2, and s 2 = [1 + fi _1 (fi r + i 7 )]/2. 

Using the mode functions Eqs. (fl"4|) and (| 1 5(1 one can calculate the correlators of the detectors for the covariance 
matrix V [l[, each splitting into two parts (( .. ) and ( .. ) v ) due to the factorized initial state. Because of symmetry, 
one has ( Q 2 A ) = ( Q B ), {P%) — ( P% )j an d ( Qa, Pb) = (Qb,Pa)- So only six two-point functions need to be 
calculated for V. 

Since q n ~ [ 7 (t — nd)/£ld] n e~~ , ( t ~ nd ^ /n\ for large t, g„ will reach its maximum amplitude (ss (n / e£ld) n / n\) around 
t — nd ~ n/^f, which makes the numerical error of the long-time behavior of V difficult to control. Fortunately 
for the late-time behavior for all d and the long-time behavior for very small or very large d, we still have good 
approximations, as we shall see below. However, before we proceed, the issue of instability should be addressed first. 

IV. INSTABILITY OF LOW-FREQUENCY HARMONIC OSCILLATORS 

Combining the equations of motion for q A A ^ and , one has 

(d 2 + 2 7 a t + n») q f\t) = ±%¥\t - d). (is) 
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where q±{t) = q A ^'(t) ± q^'fy). For t > d and when d is small, one may expand q [ ±\t — d) around t so that 



(A), 



2 7 



q¥\t)~dd t q¥>(t) + ^dhT>(t) 



or 



d 2 



Hi 
d 

27 

d 



(t) = 0( 7 d), 



JA) 



(t) = 0( 7 d)- 



(19) 

(20) 
(21) 



If we start with a small renormalized frequency fl r and a small spatial separation d < 2j/fl 2 with 7 d kept small 

so the 0(jd) terms can be neglected, then q^ will be exponentially growing since its effective frequency becomes 

imaginary (O, 2 — (2 7 /d) < 0), while oscillates without damping. A similar argument shows that q±*^ will have 
the same instability when two harmonic oscillators with small fl 2 are situated close enough to each other. 

One may wonder whether the 0(jd) terms can alter the above observations. In Appendix IA1 we perform a late- 
time analysis, which shows the same instability. The conclusion is, if Vl 2 < 2 7 /c?, all the mode functions will grow 
exponentially in time so the correlators ( Qi, Qj ) or the quantum fluctuations of the detectors diverge at late times. 
Accordingly, we define 

d ms = 2 1 /n 2 r (22) 

as the "radius of instability." For two detectors with separation d > di ns , the system is stable. For the cases with 

(7) ■ ■ • 

d = di ns , a constant solution for q+ at late times is acquired by (|20[) . while for d < di ns , the system is unstable. 
Below we restrict our discussion to the stable regime, £l 2 > 2 7 /<£ 



ZEROTH-ORDER RESULTS 



Neglecting the mutual influences, the v-part of the zeroth-order cross correlators read 



Q A {t),Q B {t) )f 



-Re- 



{ [Tl + e" 27t (fi + 2 7 e i0t sin fit)] S d 



P A (t),Q B (t) 



\(0) 



nn 2 d n + ij 

e-T* [(fl cos Qt + 7 sin M) (S d _ t + S d+t ) + (Q + » 7 ) sin fit (C d _ t - C d+t )] } 

— %-,Re i(Q + %ri) { \Q + e~ 27 * (fi - 2 7 e 4nt sin fit)] S d - 

e-^ [(ft cos fit - 7 sin fit) («S d _ t + 5 d+i ) + (0 - 17) sin fit (C d _ t - C d+t )] } 

(Q A (t),P B (t))^ 

h'y 

irn 2 d e 



-It 



sin Qt Re {-2e(-T+ in )*5 d + + S d+t + i (C d _ t - C d+t )} 



where 



(Ci[(fi + £7)2:] + Ci[-(fi + «7)x]) sin[(fi + 17)2;] - Si[(n + i-y)x] cos[(fi + i~f)x], 
(Ci[(fi + t7)x] + Ci[-(fi + «7)x]) cos[(fi + 17)2;] + Si[(fi + i~/)x] sin[(0 + 17)2;], 



(23) 

(24) 

(25) 

(26) 
(27) 



with sine-integral Si(x) = si(x) +7r/2 and cosine-integral Ci(x) [23] ■ The a-part of the zeroth-order correlators as well 



as the two-point functions (for a single inertial detector), ( Q 2 )^\ s , , . , ; 

the spatial separation d [for explicit expressions see Eq. (25) in Ref. [1| and Appendix A in Ref. [2l|]. So the d 



Qjjtj ) { °\ and ( Pf } { ° } are all independent of 



dependence of the zeroth-order degrees of entanglement E$ and E(°) are all coming from (f2"3")) - 
to the phase difference of vacuum fluctuations that the detectors experience locally. 
Note that when 

d-> d ■ = ie 1 " 7 *"^ 

where 7 e is the Euler constant and Ai = — lnf2At — j e corresponds to the time-resolution At of our detector theory 
one has ( K A (t), K B {t) )|° } -> ( K A (t) 2 ) { ° ] = ( K B (t) 2 ) { ° } , 1l = P,Q. That is, the two detectors should be seen 
as located at the same spatial point when d ~ d m in in our model, which is actually a coarse-grained effective theory. 
Let us call d m in the "merge distance." 



which are due 



(28) 
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FIG. 1: The zeroth-order results, no mutual influence is included here. 7 = 10 -5 , Q = 2.3, Ao = Ai = 20, and (a, 0) = (1.1, 4.5). 
The two plots on the left are for the zeroth-order E$ which is seen to decrease and disentangle in time. (The behavior of E*- ' 
is similar to —E$ but the amplitude of oscillation in time is smaller.) The two plots on the right are for the relative values 
of E$ at spatial separation d to the value at infinite spatial separation, as given in (I35|l . In the upper-right plot, the brighter 
color corresponds to the higher value of Ej^ xe \. 



A. Early-time entanglement dynamics inside the light cone (d < i) 



In the weak-coupling limit (7A1 <C fi), when the separation d is not too small, the effect from the mutual influences 
comes weakly and slowly, so the zeroth-order correlators dominate the early-time behavior of the detectors. The 
asymptotic expansions of sine-integral and cosine-integral functions read [22| 



Ci[(f2 + ii)x] 
Si[(f2 + ij)x] 



m / x 

"2VN 

7T X 

2R~ 



1 



sin(f2 + 17)2; 
(fi + i^)x ' 
cos(f2 + ij)x 



(Q + ij)x ' 

for fl, 7 > 0, and |(f2 + ij)x\ 1. So in the weak-coupling limit, from t — d = up to t — d ~ 0(1/7), one has 



(QA(t),Q B (t) ) { ° ] *e(t-d) 



sin fid h 

rid m l 



-yd 



1 



-27(4— d) 



(29) 
(30) 

(31) 



( P A (t),P B (t) )f pb n 2 ( Q A (t),Q B (t) >< 0) and ( P A (t),Q B (t) )f , ( Q A (t),P B (t) ) { ° ] ~ Ofr/Sl). The 9{t-d) implies 
the onset of a clear interference pattern (~ sin fld/ild) inside the light cone, as shown in Fig. [T] This is mainly due 
to the sign flipping of the sine-integral function Sid-t in (j2"3"]) - (|2l)|) around d = t when d — t changes sign. The 6(t — d) 
acts like that each detector starts to "know" the existence of the other detector when they enter the light cone of each 
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other, though the mutual influences are not considered here. In the next subsection we will see that there exists some 
interference pattern of 0(7) in £ even for d > t, where no classical signal can reach one detector from the other. 



B. Outside the light cone (d > t) 

Before the first mutual influences from one detector reaches the other, the zeroth-order results are exact. From 
and when d > t and |fi + ij\(d — t) 3> 1, one has 



<Q^(t),Q B (t))t 0) 



20 



7rfi 4 d 2 



1 + e~ 27 * (cos fit + ^ sinfiij - — — — ^cosfii — sin <>/ 



1 

n 



/ P M P (+\ \(°) ~ 27 -2 7 t sin2 
(P A (t),P B (t)) v « — e 7 -^, 

(PaW,Qb(*))1 0) = (Q A (t),P B (t))^ 



27e~ 7 ' sinffi 
TrfiHd 2 fi 



~ 7t (cos fit + ^ sin fit 



tf 2 -i 2 



(32) 

(33) 
(34) 



which makes the values of £V an d £ depend on (i and t\ that is, the dependence of the degree of entanglement on the 
spatial separation d between the two detectors varies in time t, even before they have causal contact with each other. 
In the weak-coupling limit, with the initial state ([2|) and fi 3> jAj > 7, j = 0, 1, one has 



-log 2 2c_(M)- [-log a 2c_(t,oo)] 



7 n X 



7rln2 |Af| ^ 

1 1 n=0 7 



J2 £a cosnfit + O( 7 2 A , 7 2 Ai) 



(35) 



when d > i and |fi -I- ij\(d — t) 3> 1, where 

a 7 = d" 2 {?i 2 /3 2 +a 2 (-a 2 /3 2 fi 2 +/3 4 + 4/3 2 ?ifi-n 2 fi 2 ) + |A"|(a 2 fi 2 - 

2/3 2 e- 2 ^ [?i 2 + a 2 (/3 2 - 27ifi) - \X\] } e' 27 ', 
a 7 = -4(d 2 -t 2 )-\3 2 {2a 2 ftfi + [fr 2 +a 2 (/? 2 -2?ifi)- \X\] e" 27 *} e" 



■7* 



a] - d- 2 {4hila 2 p 2 + [P 2 h 2 + a 2 (a 2 f3 2 il 2 + /3 4 - 4/3 2 ?ifi + ^ 2 fi 2 ) - |A"|(a 2 fi 2 + p 2 )] e~ 27 *} , 
6 7 = fi 3 {2h 2 na 2 p 2 + h [a 2 (a 2 /3 2 n 2 +fi 4 - 4 / 3 2 MT) + (a 2 fi 2 + (3 2 )(h 2 - \X\j\ e^ 7 ' 
+(h - a 2 fi)(^fi - (3 2 ){a 2 (3 2 + ti 2 - |) e — } , 



(36) 
(37) 
(38) 

(39) 



with X = h 2 — a 2 /3 2 . So for X ^= the relative degree of entanglement at separation d to those for the detec- 
tors at the same moment but separated at infinite distance oscillates in frequency fi and/or 2fi, depending on the 
values of a 7 . This explains the (cosfit)/(eP — t 2 ) pattern outside the light cone in the upper-right plot of Fig. [1] 
and the small oscillations before t « 7.5 in the lower-right plot of the same figure, where (a,0) = (1.1,4.5) so 
(al,al,al)/W w (1.94/d 2 , -2.89/(d 2 - t 2 ), 0.95/eP) at early times. Another example is, when (a,0) = (1.5,0.2), one 
has (<2q, aj, a^)/b 7 w (— 4.68/cZ 2 , — 0.06/(d 2 — t 2 ), 4.74/d 2 ) at early times, so the d~ 2 cos2fii pattern dominates at large 
d in the bottom-right plot of Fig. [7] For these cases, the larger the separation, the weaker the entanglement (in terms 
of the logarithmic negativity) at some moments, but the stronger the entanglement at other moments. 

The sudden switching on of interaction at t — in our model will create additional oscillation patterns outside the 
light cone. However, as shown in (|35|) . those oscillations are suppressed in the weak-coupling limit by 0(7Ao) of the 
above results. Here A = — In fi At — "f e with At corresponds to the time scale of switching on the coupling between 
the detectors and the quantum field (see Sec.IIIB in Ref.[2l[ for details). 

When f3 2 — ft 2 /a 2 or X — 0, the detectors are initially separable and 



16a 4 7r 2 fi 4 



{7rfi(7i-a 2 fi) 2 e- 27 *(l-e- 27t )+ 



2 7 Ai [2^a 2 fi(l - e~ 27t ) + /i 2 e- 27 *(l - cos2fii) + a 4 fi 2 e" 27 *(l + cos2fi£)] } + 0( 7 2 ), 



(40) 



outside the light cone, which is always positive so the detectors are always separable. When we increase the coupling 
strength 7, we find that the values of S are pushed further away from those negative values of entangled states. In 
Appendix [B] we also see that quantum entanglement is only created deep in the light cone. Therefore in our model 
we see no evidence of entanglement generation outside the light cone. 
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For \X\ but sufficiently small, the detectors are initially entangled, but after a very short-time scale 
0(e~ 7< =~( A °/ 2 )) the value of £ jumps to {h 2 ^ / (3n) 2 — (hX/4a[3) 2 , which could be positive so that the detec- 
tors become separable. In these cases quantum entanglement could revive later as E is oscillating with an amplitude 
proportional to 7A1 , while these revivals of entanglement do not last more than a few periods of the intrinsic oscillation 
in the detectors. 



C. Breakdown of the zeroth-order results 



At late times i > 7 1 , all 
read 



vanish, so 



dominate and the nonvanishing two-point correlation functions 



QA,QB) {0) \ t »^ 



Re 



Pa,Pe 



\(o) 



lt»7" 



Q 



2 v(0) I 
A I l*»7" 



Q 2 



(0) 1 
B I lt»7" 



h 
ih 



Re(in - j)S d 
ln 7 ~^ 



P 



2 \(0) 



lt»7 



-i = I P 



2 \(0) 



lt»7" 



from p3] ) -([25 |) and from Ref. [2J. 

When d — > 00, the cross correlators vanish and the uncertainty relation reads 



TO^-^det 



v(°) 



t»7" 



2 



Q 



AO) 



P 



2 \(0) 



lt>7" 



2Al -ln(l + ^ 



h 2 \ 
-r >0, 



for sufficiently large Ai [2lj], so the uncertainty relation holds perfectly. However, observing that \Sd\ 
large enough but still finite, the late-time T^ ^ can reach the lowest values: 



(QA(t)) i0) (P 2 A(t)) i0) | t »7 

h 2 



h 2 e- 2 ^ d 
Ad 2 



(o) 



lt»7" 



n 4 e -^ d 
~ wnfd 4 

n^({p 2 (t) 



(0) 



lt»7" 



(41) 
(42) 
(43) 

(44) 



(45) 

■ne~ ld for d 



(46) 



This zeroth-order result suggests that the uncertainty relation can fail if d is not large enough to make the value of 
the second line of (|4"6"]) overwhelmed by the first line [see Fig. [5]. When this happens the zeroth-order results break 
down. Therefore to describe the long-time entanglement dynamics at short distances d the higher-order corrections 
from the mutual influences must be included for consistency. 

When 7 <C 7A1 <C fi, one has a simple estimate that the late-time T*- -* becomes negative if d is smaller than about 
do w 7r/2Ai7, which is much greater than di ns found in Sec lIVl 



VI. ENTANGLEMENT AT LATE TIMES 

Since all q!f' vanish at late times in the stable regime (see Appendix A), the late-time correlators consist of q^ 



only, for example, 



Qh ' li " = I j0^&HtMQ£\t*)\^°o, (47) 



where (t,Js)\t-,oo is given by (|A12[) and q^ B ^ (t, 1c) j ^—,-00 is its complex conjugate. After some algebra, we find that 
the value of the nonvanishing correlators at late times can be written as 

( Ql ) |t-oo = ( Qb ) |t-oo = 2Re(^ 0+ +^o-), (48) 

(Qa,Qb)|^oo = 2Re(^ +-^o-), (49) 

( Pi ) l^oo = ( P 2 ) l^oo = 2Re(^ 2+ +T 2 -) , (50) 

{P A ,P B )\t^oo = 2Re(^ 2+ -^ 2 _), (51) 
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FIG. 2: The oscillating curve represents the value of T' - 1 (defined in JUJ) as a function of d. The bottom curve represents its 
lower bound (Eq.(|46}). It becomes negative when d < 616, which signifies the violation of uncertainty relation. To rectify this, 
one needs to add on the mutual influences, as shown in Fig. [3] Here 7 = 10 -4 , Q = 2.3, Aq = Ai = 25. 



where 



duj- 



2i-fcu - n 2 ± £l e 



27 „iu>d ' 



and uj m ax is the high frequency (UV) cutoff corresponding to Ai. 
In the stable regime one can write T c ± in a series form: 



(52) 



47T 



w 2 + 2*70; - U 2 — 7* 



E 



x 27 iwd 



L<j 2 + 2i~foj — fl 2 — 7 2 



J n=0 



,-_J_„iijjdn 



uj 2 + 2i-fu - SI 2 - 7 2 



(53) 



so we have 



7i i 7 — in 

— < — In 

4?r I 20 7 + in 



7; 



^ n! t 



7 



' Re e (T+«)^r[0, (7 + ifi)nd] 



— <| — (n 2 - 7 2 ) ln^S 



+ 7 



2A 1 -ln(l + ^ 



(54) 



(55) 



for large frequency cutoff w maa: , or the corresponding Ai. 

Substituting the late-time correlators (|48]) - (f5"l"j) into the covariance matrix V, we get 



( 16RaFo+Re^ 2 - - — j ^16ReJ r _ReJ r 2+ - — 

h 2 \ ( h 2 
16ReJP +ReJ r 2 + - -j J ( 16Re^b-RaF 2 - - — 



(56) 
(57) 



Numerically we found that l6ReJ-o+ReJ-2~ — (h 2 /4) and T^oo are positive definite in the cases considered in this 
paper. We then identify the late-time symplectic spectrum (c + ,c-)\ t ^oo = (4yRx^o+Rc^2^, 4yRo^oZRa^2+). So 
if lGRe^o+Re^- — {h 2 /A) is negative, then S < 0, £V > 0, and the detectors are entangled. 
In the weak-coupling limit, keeping the correlators to 0("f/d), we have 



16R&Fo+ReJ-2- 



T, 2 



?i 2 7A 1 

?rn 



n 3 



Re 



i 7 n 27 2 A X 



nd 



TT 2 d 



(i + nd) 



e lild T[0,ind] 



(58) 



11 



-2-10 



FIG. 3: Plots for E (solid curve) and T (dashed curve) at late times as a function of d, with parameters the same as those in 
Fig. [2] Two detectors are separable when E > (shaded zone). One can see that E becomes negative when d < 0.025. With 
the mutual influences included, the uncertainty relation [see Eq.© and below] now holds for all d. 



which is positive as d — > oo, but negative when d — > 0+. So {55J must cross zero at a finite "entanglement distance" 
dent > 0, where S = 0. For d < dent-, the detectors will have residual entanglement, while for d > d en t, the detectors 
are separable at late times. 

For small 7, d ent is almost independent of 7. We find that when 7A1 <C ft and Aj ^> 1, 



7r/2ft 



Ai-lnjJ- 



(59) 



will be a good estimate if d en t -C 1. Here d en t is still much larger than the "merge distance" d m i n in (|28|) . For 
example, as shown in Fig. [3l when 7 = 0.0001, ft = 2.3, Ai = 25, one has d ent ~ 0.025, which is quite a bit greater 
than the "radius of instability" 27/ft 2 w 3.8 x 10 -5 , and much greater than the merge distance d m in ~ 9 x 10 -12 . 

A corollary follows. If the initial state of the two detectors with d < d ent is separable, then the residual entanglement 
implies that there is an entanglement creation during the evolution. In contrast, if the initial state of the two detectors 
with d > d en t is entangled, then the late-time separability implies that they disentangled in a finite time. Examples 
will be given in the next section. 

Note that the ill behavior of has been cured by mutual influences. The uncertainty function ([57]) is positive 
for all d at late times. 

Note also that, while the corrections from the mutual influences to ( Q 2 A ) \t->oo and ( P\ ) |t-«x> are Oi^y/d), the 
mutual influences have been included in the leading order approximation for the cross correlators. Indeed, in (|53D . 
even as low as n = 1, we have had 



Qa, q 



B / |t^oo 



Qa,Ql 



AO) 



2h"i 47 

7T d 



dco- 



[(fi 2 . — uJ 2 ) cosuid — 27wsinu;<i] 
[{to 2 - ft 2 .) 2 + 4 7 2 u; 2 ] 2 



(60) 



However, this is slightly different from the approximation with the first-order mutual influences included. Writing the 
n = and n = 1 terms in Eq. (fl4|) as 



,(+) 



o (+) W +) 



(61) 



then the approximated cross correlator with the first-order mutual influences included is the u integration of Re 
[(?£Lo + <?K=i)(<4tn=o + QbI=i)}> but in © only Re [q^ n=0 qff n=0 + Q^LJeLi + ^Li^Lol contribute, 
though there are 0(7°) terms in <ZAn=lSs n=i- ^ n e latter is small for fid 3> 1, and will be canceled by the mutual 
influences of higher-orders. 



VII. ENTANGLEMENT DYNAMICS IN WEAK-COUPLING LIMIT 
A. Disentanglement at very large distance 



Suppose the two detectors are separated far enough (d 3> ft) so that the cross correlations and the mutual influences 
can be safely ignored. Then in the weak-coupling limit (ft ^> 7A1) the zeroth-order results for the v-part of the self 
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correlators dominate, so that [l[ 

{<&). = < Q% >v - i (1 " ^) - (62) 

( Pi ) v = ( Pi K « ^ (1 - e"^) + ^ 7 A 1; (63) 

and ( Qa, Pa ) v = ( Qb, -Pb ) v ~ 0(7), while the v-part of the cross correlators are vanishingly small. This is exactly 
the case we have considered in Sec. IV A 2 of Ref. [l|, where we found 

with Z$ > 0, Z% — Z4 > and Z2 > [Zg, Z4 and Z2 are parameters depending on a and f3, defined in Eqs.(37), (38) 
and (41) of Ref. [l|, respectively] Accordingly the detectors always disentangle in a finite time. There are two kinds 
of behaviors that £ could have. For Z4 > 0, the disentanglement time is a function of Z4, Z% and 7, 



while for Z4 < 0, the disentanglement time is much longer, 



dB< 2 7 Z 2 + - 4a 2 /3 2 Z 4 V ' 



and depends on Ai. 



B. Disentanglement at large distance 



When d is large (so l/fid is small) but not too large to make all the mutual influences negligible, while the zeroth- 
order results for the v-part of the self-correlators (f6"2"]l and (j6"3")) are still good, the first-order correction [n = 1 terms in 
p4p] to the cross correlators ( Qa, Qb ) can be of the same order of ( Qa, Qb )^ (a similar observation on the late- 
time correlators has been mentioned in the end of Sec. IVI[) . Including the first-order correction, for d > 0(l/\/ 7 f2), 
we have a simple expression, 



Qa, Qi 



{Qa, 
0(t- 



d) 



h sin fid 



d) 



n rid 



2Q, nd 

e^ d 1 (t - d)e- 2 ^*- d) . 



-1 + e - 2 T(*- d ) (1 + 2(t - d) 7 ) + 0(7/n) 



(67) 



\(0) 



111 




and ( P4, Ps ) v ~ " ( Qa,Qb ) v with other two-point functions ( .. ) v being 0(7) for all t. Here ( Qa,Qb 
the weak-coupling limit has been shown in (|3ip . The above approximation is good over the time interval from t 
up to e - 2 ^*- d ) > 0(7/0), namely, before t - d ~ 0(~7 _1 lnfr/fi)). 

Still, in this first-order approximation, ( Qa,Qb ) v and ( Pa,Pb ) v are the only correlators depending on the 
separation d. Inserting those approximated expressions for the correlators into the definition of S or Ej^f, we find 
that the interference pattern in d for the relative values of £ or Ej^f at early times (Fig. [T]) can last through the 
disentanglement process to make the disentanglement time tdE longer or shorter than those at d — > 00, though the 
contrast decays noticeably compared with those at early times. Two examples are shown in Fig. |4] For Z4 > 0, the 
disentanglement time is about 



tdE^ 



b dE> 





z 6 ( 


V°) 

b dE> ' 




e 7d sin fid 




Z s d( 


L-e- 2 <i>) 


+ z 6 


1 - 


^(t { d°L-d)\ 


e< d sin fid 



(68) 



where Z 6 ee (h 2 - a 2 /3 2 )(a 2 fi 2 - /? 2 ) [Fig. H (left)]. In this case the disentanglement time can be short compared 
to the time scale 0(71/7), n £ N when the higher-order corrections q n from mutual influences reach their maximum 
values (see Sec. IIII[) . So in the weak-coupling limit the above estimate could be good from large d all the way down 
to fid ~ O(l) but still much greater than fid en t- If this is true, the difference of disentanglement times for different 
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FIG. 4: The plot of E as a function of d and t, up to the first-order correction. E is negative in the dark region and positive 
in the bright region. For a fixed d, the disentanglement time tdE is at the border of the lowest dark region or the earliest time 
that the detectors becomes separable. The interference pattern in Fig. [T]for E at early times signifies that the disentanglement 
time tdE is longer or shorter than those at d — > oo [Eqs. (|65| and (|66|l 'l. The gridded profile in the left plot shows that after 
tdE there could be some short-time revivals of entanglement. Here the parameters are the same as those in Fig. [1] except 
(a, /3) = (1.5,0.2) in the left plot and (1.1,4.5) in the right (cf. Fig. 3 in Ref. [U). 

spatial separations can be significant at small d. For example, for (a, (3) — (1.5,0.2) with other parameters the same 
as those in Fig. [5J the disentanglement time at d 4.4934/fi (where smftd/ild is the global minimum) is over 1.6 
times longer than those for d ~ 7.7253/il (where the first peak of sin ild/ftd is located). 

For Z4 < 0, the correction of sinfM is below the precision of estimated in (jSB")) . Here we just show the 

numerical result up to the first-order mutual influences in Fig. U (right), which shows that the interference pattern in 
d is suppressed but still nonvanishing for large disentanglement times. 



C. Entanglement generation at very short distance 



When fid ~ O(e), "f/Cl ~ 0(e 2 ), and e <C 1, one can perform a dimensional reduction on the third derivatives in 
(TT^l) . namely, 



to obtain, up to 0(e 5 ), 



qf + 2 1± q { £> + n 2 ± q { £ 



ll r + d 

IT id 



0, 



where j = A, B, q { ± ] = <?^ +) ± q { + ] and 

_ 7d2 (ni + 22) 



6 (l+ 7 d) 2 
2 7 1 d 2 (0 



'2 fY)2 _ 27 \ 
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1-jd 6(1- 7 d) 2 ' 
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FIG. 5: (Upper left) The solid curve and the long-dashed curve represent the values of E and T, respectively, while the dotted 
line is for the value of E with all { .. ) v set to zero. The detectors are separable initially (the parameters here are the same as 
those in Figs. [2] and [3] except d = 0.01 and {a, 13) = (1, 1)). Quantum entanglement has been generated after t ~ 0.15, and E 
oscillates in frequency Q+ at early times. Around the time scale t ~ 1/7+ (~ 5000 here) (upper right), E has an oscillation 
with long period 7r/(fi_ — fi+) (~ 361.4). E appears to be settling down at a value (about —0.046 here) depending in this case 
only on the value a in the initial data. However, in a much longer time scale t ~ I/7- (~ 1.13 x 10 8 ) (lower left), one sees that 
the value of |E| is actually decaying exponentially to the late-time value (w —6.8 x 10~ 7 ) consistent with the results in Fig. [3] 
with d = 0.01 and independent of the initial data of the detectors. 



Here 7-/7+ is of 0(e 2 ). Note that q>_ and the decay modes in have subradiant behavior, while q 3 + and the decay 

modes in q^&' are superradiant. For small d, the time scale ^Z 1 3> 7 _1 > 7+ 1 ~ 1/27, an d 7Z 1 goes to infinity as 
<2->0. 

The solutions for (|70| and (f7Tj) with suitable initial conditions are 



,0) 



±3 



U) _ 



a (+) ± n (+) 
1a ± 1b 



A 



± c -iQ±t] 



n+ V 



1.1/2 



±e 



s, e 



sfcid/2 



[(M± 



Mjf ) e~ lwt + e- 7±t (M 



± if2+t 



)] 



(75) 
(76) 



where = [l-fi± 1 (fi r +i7 ± )]/2, s± = [l+ri± 1 (ri r +i-f±)}/2, Mf = (-uj-i-/ ± +n ± )- 1 , and Mf = (-uj-i-/±-tt ± )- 1 . 
Actually these solutions are the zeroth-order results with 7 and Q replaced by 7± and Q±. So we can easily reach the 
simple expressions 



167T7 + 

A 2 + 

167T7+ 



Qa, Qb 



Qa, Qb 
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f + ( Qa > v 0) 



n->n + 
7^7+ 
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167T7_ 

A 2 _ 

167T7- 



Qa )f -<QA,Q B )f 



v(o) 



v(0) 



o->n_ 

7^7- 



(77) 



(78) 



and so on. Here 

({Qa,Pa) v = d t 



3 above and in Eqs.(A9) and (A10) of Ref.[21| 



)^ are those expressions given in 

\ ) v /2.) The prefactors A? t /167r7± are put there because in our definitions for the zeroth- 
order results the overall factor A 2 , has been expressed in terms of 8777, but now j± ^ X±/8tt. 

In Fig. [5] we demonstrate an example in which the two detectors are separable in the beginning but get entangled 
at late times. There are three stages in their history of evolution: 
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1. At a very early time (t ~ 0.15) quantum entanglement has been generated. This entanglement generation is 
dominated by the mutual influences sourced by the initial information in the detectors and mediated by the field. 
(For more early-time analysis, see Appendix 151) 

2. Then around the time scale t ~ I/7+, the contribution from vacuum fluctuations of the field (( .. ) ) takes over 
so that £ becomes quasisteady and appears to settle down at a value depending on part of the initial data of the 



detectors. More explicitly, at this stage q+ 
their initial values, so 



O) 



fi = A, B, +, — have been in their late-time values but q_ are still about 



sin ^ -2 7 d , , 2 

fid + n 



sin fid _ 27d _ 2h_ 8A17 
fid 6 a 2 fl nfl 



(79) 



in the weak-coupling and short distance approximation 7 <C dfl 2 <C fl. Here £ depends on a only. The parameter 
P in initial state ([2]) is always associated with qr) in ( .. ) so it becomes negligible at this stage [cf. Eq.(25) in [l[]. 
Note that £| t/sjl / 7+ can be positive for small d only when a is at the neighborhood of y/h/fl. 

3. The remaining initial data persist until a much longer time scale t ~ l/-f_ when £ approaches a value consistent 



with the late-time results given in Sec. IVI1 which are contributed purely by the vacuum fluctuations of the field and 
independent of any initial data in the detectors. In this example the detectors have residual entanglement, though 
small compared to those in stage 2. 

The above behaviors in stages 2 and 3 cannot be obtained by including only the first-order correction from the 
mutual influences. Thus in this example we conclude that the mutual influences of the detectors at very short distance 
generate a transient entanglement between them in midsession, while vacuum fluctuations of the field with the mutual 
influences included give the residual entanglement of the detectors at late times. 

For the detectors initially entangled, only the early-time behavior looks different from the above descriptions. Their 
entanglement dynamics are similar to the above in the second and the third stages. 



VIII. DISCUSSION 



A. Physics represented by length scales 

The physical behavior of the system we studied may be characterized by the following length scales: 

Merge distance d m i n in Eg. (I28p .-Two detectors separated at a distance less than d m i n would be viewed as those 
located at the same spatial point; 

Radius of instability d ins in Eq. (I22[) .-For any two detectors at a distance less than di ns , their mode functions will 
grow exponentially in time so the quantum fluctuations of the detector diverge at late times; 

Entanglement distance d en t in i?o. (|59p .-Two detectors at a distance less than d en t will be entangled at late times, 
otherwise separable; 

And do defined in Sec. IV Cl -For d < do the zeroth-order results breakdown. A stable theory should have d en t and 
dmin greater than d ins . 



B. Direct interaction and effective interaction 



In a closed bipartite system a direct interaction between the two parties, no matter how weak it is, will generate 
entanglement at late times. However, as we showed above, an effective interaction between the two detectors mediated 
by quantum fields will not generate residual entanglement (though creating transient entanglement is possible) if the 
two detectors are separated far enough, where the strength of the effective interactions is weak but not vanishing. 



C. Comparison with 2HO QBM results 



When d — > d m i n with large enough fl, our model will reduce to a 2HO QBM model with real renormalized natural 
frequencies for the two harmonic oscillators. Paz and Roncaglia [ll| have studied the entanglement dynamics of this 
2HO QBM model and found that, at zero temperature, for both oscillators with the same natural frequency, there 
exists residual entanglement at late times in some cases and infinite sequences of sudden death and revival in other 
cases. In the latter case the averaged asymptotic value of negativity is still positive and so the detectors are "entangled 
on average." 
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While our results show that the late-time behavior of the detectors is independent of the initial state of the detectors, 
the asymptotic value of the negativity at late times in [l5T ] does depend on the initial data in the detectors (their 
initial squeezing factor). This is because in [l5| the two oscillators are located exactly at the same point, namely, 
d = 0, so 7_ = and the initial data carried by q_ persists forever. Since in our cases d is not zero, the "late" time 



151 ] actually corresponds to the time interval with (I/74.) <Ci<C(l/7-)in our cases, which is not quite late for 



our detectors. 



D. Where is the spatial dependence of entanglement coming from? 



Two factors are responsible for the spatial dependence of entanglement. The hrst one is the phase difference of 
vacuum fluctuations that the two detectors experience. This is mainly responsible for the entanglement outside the 
light cone in all coupling strengths and those inside the light cone with sufficiently large separation in the weak- 
coupling limit, such as the cases in Sec|Vj The second factor is the interference of retarded mutual influences, which 
are generated by backreaction from the detectors to the field. It is important in the cases with small separation 
between the detectors, such as those in Sec lVII Cl 



E. Non-Markovian behavior and strong coupling 

In our prior work [2l| and [l[ , the non-Markovian behavior arises mainly from the vacuum fluctuations experienced 
by the detectors, and the essential temporal nonlocality in the autocorrelation of the field at zero temperature 
manifests fully in the strong-coupling regime. Nevertheless, in Sec. IVII Cl one can see that, even in the weak-coupling 
limit, once the spatial separation is small enough and the evolution time is long enough, the mutual influences 
will create some non-Markovian behavior very different from those results obtained from perturbation theory with 
higher-order mutual influences on the mode functions neglected. 
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APPENDIX A: LATE-TIME ANALYSIS ON MODE FUNCTIONS 



'1+ 



W (t) = 5>e«", (Al) 
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Let 



Equation (TTg}) gives 

c 3 [- K 3 2 + 2 ^Kj + ^r] e^* = y Cj-j* 1 '*"*' (A2) 
3 3' 

At late times, one is allowed to perform the Fourier transformation on both sides with t integrations over (—00,00) 
to obtain 

- K 2 + 2i 1 K j + Q 2 r = ^e-^ d . (A3) 
d 

There are infinitely many solutions for Kj in the complex K plane, so one needs infinitely many initial conditions to 
fix the factors Cj. Our q + chosen as a free oscillator at the initial moment and unaffected by its own history until 
t = d in principle can be specified by a set of c^'s. Suppose this is true. Writing Kj = Xj + iyj, the real and imaginary 
parts of (|A3|) then read 

(y - 7 ) 2 - x 2 + Q 2 = ^-e yd cosxd, (A4) 

x(y-l) = -,e vd sinxd. (A5) 
a 
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FIG. 6: (Upper left) The solutions to ()A3|) for the complex frequency Kj — Xj + iyj of [defined in (|A1|) ] are located 
at the intersections of the dashed and solid curves, which represent the solutions to Eqs. (|A4[l and (|A5[) . respectively. Here 
7 = 0.25, Q — 0.9000202, d = 1. (Upper right) The same case, but here "+" denotes complex solutions and "x" denotes purely 
imaginary solutions for K. There are two purely imaginary solutions for K in this case. (Lower left) There are three purely 
imaginary solutions for K when 7 = 0.25, Q — 0.8, d — 1. (Lower right) Solutions for K when 7 = 0.25, = 0.3, d= 1. There 
is only one purely imaginary solution, which is located in the lower half of the complex K plane. 



The solutions for them are shown in Fig. \E\ The left-hand side of (|A4|) is a saddle surface over the xy space, while 
the right-hand side of (|A4[) is exponentially growing in the +y direction and oscillating in the x direction. For (|A5[) . 
the situation is similar. From Fig. [6l one can see that there is no complex solution for K with nonvanishing real part 
and negative imaginary part (i^0 and y < 0). The solutions for K with its imaginary part negative must be purely 
imaginary. Indeed, from (|A5|) and Fig. [5] (upper right), one sees that when x 7^ 0, if y < 0, then (y — 7) < —7, but 
— 0.21727 < j e v d (sin xd)/(xd) < 7, so there is no solution of (|A5|) with y < and x ^ 0. 

When > 27/ei, one finds that all solutions for K in (|A3[) are located in the upper half of the complex K plane, 
i.e., all yj > 0, which means that all modes in (|A1[) decay at late times. 

When = 27/d, there exists a solution if = 0, with other solutions on the upper half K plane. This implies that 

q + becomes a constant at late times. 

When fi^ < 27/d, there must exist one and only one solution for K with negative y, which corresponds to the 
unstable growing mode. This is consistent with our observation in Sec lIVI 

Therefore, we conclude that is stable and decays at late times only for > 27/d. 



As for qi^ , from 



T|) it seems that q 



(A) 



that q^ decays at late times for all cases. Thus, by symmetry, all q^ 
fl 2 r > 27/d. 



would oscillate at late times. However, similar analysis gives the conclusion 

decay at late times in the stable regime 
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Now we turn to q^ B - Equation (fl3l) implies that 

(a i 2 + 2 7 a t + a?)%W(t,k) = qB +\t- 2 d,k) + 

(- W 2 -2 7W + f> 2 )e^ d / 2 



at late times. Again, let 



A e" 



<4 +) (i,k)=]T<e 



iKlt 



2 7 



iujd—ik\d/2 



(A6) 



(A7) 



then one has 

£4 



2»7«/ + U 2 



\ e 



\ _ — iujt 
Ar 



(~uj 2 - 2iju; + fi 2 ) e ifclrf / 2 + ^e Md - jfelrf / 2 



After a Fourier transformation, for Ki ^ —lu, the above equation becomes 



Kl) + 2 i7 i^ + f! 2 



which is the square of Eq. (|A3[) for qY^ 1 , or the square of the counterpart for q^ . So these modes decay at late 



?1) e -2H<;..i 
d 



(A) 



(A8) 



(A9) 



times for Vt 2 > 2 7 /d as q^ and q 1 ^' do. On the other hand, if, say, K 1 ^ = —to, one has 



—Hud 



Ac 



(-c^ 2 - 2z 7 ^ + n 2 r ) e^ d ' 2 + 



This equation will not hold unless 



„ _ A [{-uj 2 - 2ijuj + fl 2 r ) e lk ^/ 2 + 
Ck [-lu 2 + 2i~fLU + n 2 r } 2 ~ (£ff e -^d 

Therefore, for il 2 > 2j/d, the only mode which survives at late times will be e~ lut , and 

(+) / / i \ I —iu)t 

TB (^ k )|t»l/7 = C k e 



(A10) 



(All) 



(A12) 



This is nothing but the sum of the e -M*-" d ) p ar t i n Eq. (fT4j) with t — > oo so summing from n — to oo. Thus, (|A12[) 
with (jAllj) has included the mutual influences to all orders. The above analysis also indicates that the e _ T( t_n(i ) part 
in (fT4|) really decays at late times for f2 2 > 2 7 /d. 



APPENDIX B: EARLY-TIME BEHAVIORS IN WEAK-COUPLING LIMIT 



In the weak-coupling limit, the cross correlators ( 1Za,71 b ) with 1Z. 7Z' — Q,P are small until one detector enters 
the other's light cone. From this observation one might conclude that the cross correlations between the two detectors 
are mainly generated by the mutual influences sourced by the quantum state of the detectors and mediated by the 
field. This is not always true. 

As shown in Sec IV A[ the interference pattern inside the light cone has been there in the zeroth-order results, where 
the mutual interferences on the mode functions are not included. A comparison of the first-order results in the upper 
plots in Fig. [7] and those of the zeroth-order in Fig. [1] shows that the corrections to entanglement dynamics from 
mutual influences at early times are pretty small in that case. Actually the early-time dynamics of entanglement in 
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FIG. 7: The early-time evolution of £V with the first-order mutual influence included for different initial states of the detectors 
with 1/8 < d < 15. (Upper row) All parameters are the same as those in Fig. [T]where (a,/3) = (1.1,4.5) and the initial state 
of the detectors is entangled. Compared with Fig. [T] one can see that the distortion of the interference pattern due to the 
mutual influences is tiny. (Lower row) (a,/3) = (1.5,0.2), the detectors also initially entangled. The distortion by the mutual 
influences is also tiny. As indicated by Eq. (|35p . the complicated structure of Ej^ Te \ outside the light cone is reducing to simple 
oscillations as time goes larger. 
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FIG. 8: The early-time evolution of £V for an initially separable detector pair with the first-order mutual influence included. 
The parameters are the same as those in Fig. [1] except (a, (3) = (1,1) and 1/15 < d < 15 here. From the left plot, one 
finds that quantum entanglement is created at small d due to the mutual influences. In the right plot, one can see that 
there is no clear interference pattern in d similar to those in Fig. [T]for Ej^ Ye \ inside the light cone. Note that here -EVrci = 
— log 2 2c_ (t, d) + log 2 2c- (t, oo) instead of Etf(t, d) — Ej\f(t, oo). The detectors with smaller separation d always get a greater 
value of — log 2 2c_ . While there are small oscillations outside the light cone, the smaller separation d always associates the 
greater value of .E/Vrel- The amplitude of the small oscillation is the same order of 7 2 Aq and j 2 A\. 



20 



both examples in Fig. [7]are dominated by the zeroth-order results, thus by the phase difference of vacuum fluctuations 
in ( IZatTZ-'b )v rather than mutual influences. One can see this explicitly by inserting the mode functions in the 
weak-coupling limit with the first-order correction from the mutual influences into Eq.(25) in Ref. and write 



E(t) m E + <r[ 0) t + a^>e + 6{t - d) a\'> {t - d) + a^> {t - df + 0( 7 a ) 



JO) +2 



J 1 )/ 



(Bl) 



at early times when 0(e~ 7e ~ (Ao/2) /^) < t < 0(l/ 7 Ai), z = 0, 1. Here E , cti, and cr 2 depend on a, /3 and of 0(7°), 
0(7), and 0(7 2 ), respectively Then it is easy to verify that mutual influences are negligible in the dominating 
term after 8(t — d) for the initial states with the value of 1 not in the vicinity of h jo? or a 2 Q 2 . 

In contrast, if the initial state is nearly separable {(3 2 « h 2 /a 2 ), mutual influences will be important in the 
detectors' early-time behavior. In this case, dropping all terms with small oscillations in time, the factors in (|Bip are 
approximately 



So 



JO) 



Ji) 



JO) 



r 



[n 2 7 Ai+a 4 fi 2 7(2A + Ai)] ; 



in 2 



7 



2 h 2 (h 2 



2irfl 3 



4f2 4 a 4 d 2 
2ft 2 7 A - 



& 7A1 



(?i-a 2 ft) , (B2) 



with <Ti negligible. So E evolves as the following. In a very short-time scale 0(e~ 7e ~( A °/ 2 )/il) after the interaction is 
switched on, E jumps from its initial value (ss 0) to a value of the same order of Eo, which is positive and determined 
by the numbers A and Ai corresponding to the cutoffs of this model (the difference from the exact value is due to 
the oscillating terms dropped). For a 2 ^ h/fl so Qa and Qb are each in a squeezed state initially the detectors keep 



separable at t < d since af ^ and a^'' are positive definite. But a^ 1 is negative and proportional to l/d 2 , thus after 
entering the light cone of the other detector, if the separation d is sufficiently small, or 



JO) 



(i) 



d < di 



h + a 2 n 



2 n 



(B3) 



can overwhelm and alter the evolution of E from concave up to concave down in time. If this happens, the 



quantity E could become negative after a finite "entanglement time" 



jo) 



+ 0-. 



(i) 



jo) 



2ai 1] d + 



Or 



- 2ai 1] d 



JO) 



J 1 ) 



Ji) 



rf 2 



(B4) 



This explains the entanglement generation at small d in Fig. [51 [Note that the above prediction could fail if t ent > 
0(l/ 7 Ai), i = 0, 1, and even for t ent < 0(l/ 7 Aj) the above estimate on t ent could have an error as large as 0(2n/Q) 



due to the dropped oscillating terms.] The first-order corrections to ( .. ) contribute the cos 2 fid part of a^' = 

(cos 2 fid + sin 2 fid) , so for those cases with separations small enough such that sin 2 fid -C cos 2 fid the early-time 
entanglement creations are mainly due to mutual influences of the detectors, which is causal. 

d\ in (|B3[) can serve as an estimate for the maximum distance that transient entanglement can be generated from 
a initially separable state in the weak-coupling limit, while for the detectors with the spatial separation between d\ 
and dent the transient entanglement generated at early times will disappear at late times. 
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